Pseudo-spin canting transition in bilayer quantum Hall ferromagnets: a self-charging 

capacitor 



O 

o 

(N 

U 

Oh- 
< 



D 

a 

•4-3 

a 

I 

C 

o 
o 



> 

00 

o 



o3 

s 

I 

c 

O 
o 



Two of many striking and experimentally testable con- 
sequences are (i) the universal nonlinear behavior of the 
interlayer charge imbalance q(V, B») with gate voltage V 
and in-plane field Bu 



q(V,B 
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for fcii > b c (g), leading to a spontaneous interlayer charge 
imbalance q, even in the absence of an applied interlayer 
voltage, V — ► 0, i.e., the QH bilayer is a self-charging 
capacitor, and 
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For sufficiently strong in-plane magnetic field a vt = 1 bilayer quantum Hall pseudo-ferromagnet 
is expected to exhibit a soliton lattice. For sufficiently close layers and large in-plane field, we 
predict this incommensurate "planar" phase Pj to undergo a reentrant pseudo-spin canting tran- 
sition to an incommensurate state Cj, with a finite out-of-plane pseudo-magnetization component, 
corresponding to an interlayer charge imbalance in regions between solitons. At T > the tran- 
sition is in the 2d compressible Ising universality class, and at T = 0, the quantum transition is 
in heretofore unexplored universality class. The striking experimental signatures are the universal 
nonlinear charge- voltage and in-plane field relations, and the divergence of the differential bilayer 
capacitance at the transition, resulting in a bilayer capacitor that spontaneously charges itself, even 
in the absence of an applied interlayer voltage. 

PACS: 73.20.Dx, ll.15.-q, 14.80.Hv, 73.20.Mf 

There now exists considerable experimental [[l] Q and 
theoretical (|-(|] evidence for the bilayer phase coherent 
quantum Hall (QH) state at a total electron filling frac- 
tion vt — 1, which is driven by Coulomb exchange in- 
teractions and is expected to survive even in the limit 
of vanishing interlayer tunneling, A. A number of com- 
plementary pictures of the state, include superfluidity of 
Chern-Simons composite boson j7j, excitonic superfluid, 
and pseudo-ferromagnet 

In the latter picture, that we employ here, the z com- 
ponent m z = (ni — n 2 )27r^ 2 (I = y/hcfeB is the mag- 
netic length) of the pseudo-spin magnetization unit vec- 
tor rh = + m z z, is the normalized interlayer imbal- 
ance in electron layer densities 711,2, while the azimuthal 
angle <j) = 4>i ~ 4>2 of mi is the difference between the 
electron phases <p\^ hi the two layers. Interlayer charg- 
ing energy, 2ir£ 2 e c , explicitly breaks SU(2) pseudo-spin 
symmetry down to L7(l) x Z 2 , forcing m to lie in the 
easy-xy-plane, defined by m z = 0. The interlayer tun- 
neling energy, acts like a pseudo-magnetic field directed 
along x, further explicitly breaking U(l) symmetry. 

In this Letter, we predict and explore a striking, reen- 
trant quantum phase transition from the charge-balanced 
"planar" (P) state into an interlayer charge-imbalanced 
"canted" (C) state. This novel state is characterized 
by the development of m z ^ and corresponds to the 
pseudo-magnetization rh spontaneously canting out of 
the easy-xy-plane. The transition is controlled by di- 
mensionless parameters g = 27r^ 2 e c /A and 6|| = B^/Bci, 
with a T = phase boundary b c (g) illustrated in Fig|l], 
and Bqj a critical in-plane field for the commensurate- 
incommensurate (CI) transition. j|Jl^] For a range of pa- 
rameters, we find that the transition can be continuous 
and for T — is in a heretofore unexplored 2 + 1-d quan- 
tum compressible Ising universality class. At finite T this 
transition is the 2d compressible Ising universality class. 
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FIG. 1. Phase diagram illustrating a transition between the 
"planar incommensurate" Pi and "canted incommensurate" 
Ci phases of vt = 1 QH bilayer, with insets showing cartoons 
of the two phases. 

(ii) the concomitant divergence of the interlayer differ- 
ential capacitance C = dq/dV 
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C(V,B\\) = m-b c \-iC(V\b\\-b c \-*) 
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at the quantum (T = 0) P-C transition, with exponents 
(3 = 1/2, x = 3/2, 7 = x — P = 1, and the dimensionless 
scaling functions g(a;) x — x' 3 /*, C'(ic) ^ x" 1 ^ that 
are universal. Both relations should be readily testable 
in a balanced gate geometry with the QH bilayer acting 
as a novel dielectric medium. 

We also expect a universal suppression of the quantum 
Hall gap below the P-C transition, with 6Aqh(V, Bu) oc 
—q(V, -B||) 2 , accessible through activated behavior of the 
longitudinal resistivity. Similar predictions also hold for 
the finite T transition, with b\\, b c (g) replaced by T, 
T c {b\\,g), and (for an incompressible soliton lattice) with 
exact critical exponents of the 2d Ising model. 

The driving force for the P-C transition is the competi- 
tion between the charging (hard axis anisotropy) energy 



density £ cha rge 



\s c m z , minimized by m z = 0, and the 



exchange energy density ^exchange = \ps\V<f>\ 2 - Because 
p s » — nil) ( a * least in mean-field theory), the latter 
can be lowered for any planar twisted state by increasing 
m z , via pseudo-spin canting out of the xy-plane, cor- 
responding to the development of the interlayer charge 
imbalance. In the simplest A — > limit, for a uniformly 
twisted, 4> = Q x staggered current carrying state, this 
canting instability will clearly take place for Q > \J e c / 
corresponding to an estimate of the "depairing" stag- 
gered (spin) current J c s = fp s V(/) = f y/p s e c ±. |t| In 
contrast, for a finite tunneling energy A there exist equi- 
librium twisted states, induced by an applied in-plane 
field B|| = Buy. Although the general driving mecha- 
nism is the same, as we discuss below, the details depend 
on whether the canting phase transition takes place out 
of the planar commensurate Pc or planar incommensu- 
rate Pi (soliton) state. || 

We begin our analysis by studying the classical ener- 
getics of the quantum Hall bilayer in the presence of an 
in-plane (physical) magnetic field B|| = B\\y. The appro- 
priate Hamiltonian, in the absense of quenched disorder, 
and for a fixed uniform value of m z is given by 



A 

2^P 



cos 



-Qx] 



(3) 



Aq(1 — ml) 1 / 2 is the interlayer tunnel- 
I Q is the in-plane "magnetic wavevector" 



where A = 

ing energy, _ 

Q = 2-kB\\<1/4>q associated with B\\, and A = yj p s 2irl 2 /A 
is a screening length for spatial deformations of 4>(x). In 
the commensurate state, |||l0| stable for Bu < Bci = 
2</>o/(7T 2 (iA), the pseudo-spin rh twists uniformly about 
z-axis with <p(x) — Qx. The corresponding m z depen- 
dent energy density is given by 



£c[m z 
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where the first term is the planar (i.e., m z — 0) exchange 
energy, the second consists of the interlayer charging en- 
ergy, with b = |-^/27r£ 2 £ c /A = \g 1 l 2 and the reduc- 
tion of the exchange energy due to canting, and third, the 
tunneling energy. jll| Simple analysis of £c[m z ] predicts 
a continuous canting transition from m z = planar to 
m z ^= canted state, at 6m = bi°\g) — jy/g + 1. Given 
that it is Z 2 symmetry that is being broken, at T = 
(T 7^ 0) we generically expect such quantum (classical) 
Pc — Cc transition to be in the well-studied 3d (2d) Ising 
universality class, and in the latter case characterized by 
exactly known exponents. In this scenario, this Pc — Cc 
canting transition will be followed by the Cc — Cj tran- 
sition in a CI universality class, modified by long-range 
dipolar interactions p3 and delayed by m z ^ to higher 
critical in-plane field Bci = B CI (1 — to 2 ) -1 / 4 . 

However, because currently studied bilayer devices 
are characterized by large <?, and therefore have 



b^ig) > 1, above results suggest that the system will 
first undergo a commensurate (planar)-incommensurate 
(planar) (Pc — Pi) transition at 6|| = 1 into a state char- 
acterized by a periodic array of Id y-directed solitons. 
flfTo! Consequently, we need to analyze the Pj — Cj tran- 
sition from within such incommensurate soliton state. 

Although there has not been any direct experimental 
evidence, the CI transition picture was successfully used 
H to interpret a precipitous drop in the QH gap, upon 
application of sufficiently strong B\\ Q . In addition to the 
intrinsic interest in the Pi — Ci transition, the more eas- 
ily detected canied-incommensurate (Cj) state will also 
facilitate the detection and study of the soliton lattice. 

The properties of Id incommensurate soliton state 
are well documented in the literature. For b\\ > 1 
(Q > Qci)i the soliton chemical potential is negative and 
y-directed solitons enter the system at a density n s (b\\) 
that continuously increases with bu. In this state, the 
phase (f> s (x) rotates uniformly between solitons, but slips 
by 2ir over a soliton width A, with m±(x) failing to follow 
the winding rate Q imposed by Bu . Generically, solitons 
reduce the energy density relative to that of the commen- 
surate state with £i[m z ) = £c[m z ] - S£ so n tons [m z ]. 

Although, in the absence of fluctuations, exact expres- 
sions for 4> s (x) and the energy £/ (expressible in terms 
of elliptic integrals) are available, |1(J sufficiently close to 
the CI transition thermal fluctuations |T^ | always quali- 
tatively modify these T = predictions. Consequently, 
there are at least three regimes above the CI transi- 
tion. Asymptotically close to the transition, where the 
soliton array is sufficiently dilute, such that n s (bu) < 
n T , with n T determined by n T X = |e" 1 / (2 " !rA) w 
l/(21ogp s /T), thermal fluctuations-induced steric inter- 



action, T 2 X/(8p s \xi 
nentially weak T = 



— Xj\ ) dominates over the expo- 
interaction £iL e -\ x i- x 3 Such 
enhanced soliton repulsion leads to significantly slower 



(4) (than T = 0, n B (bu) 



| log (fey — 1)| prediction), 
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power-law increase in soliton density, n s (b\\) ~ xj\b\\ 
ll 1 / 2 , and energy density given by 



5£ 



(T) r 
solitons I- 
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Ignoring for simplicity a possible intermediate T = 
dilute regime, defined by ut ^ ^s(^li) <C A -1 , f° r suffi- 
ciently high soliton density <C n s(bu) < A -1 , soliton 
interaction crosses over to contact interaction, with en- 
ergy scale p s . In this dense regime n s (b\\) ~ A _1 |fe|| — 1| 
and soliton energy density is given by 



S£ 



(dense) r 
solitons I 



2tt£ 2 



(l-m^T- (6) 



Finally in the super-dense limit, n s (bn) 3> A , </> s (^) 
only shows periodic modulation with vanishing ampli- 
tude and the incommensurate energy density reduces to 



^(supcr-dcnsc)^ 



-ra, 



128^ 2 b 2 



(7) 



with exchange and tunneling energy only leading to a 
vanishing, m z -independent correction to the charging en- 
ergy. 

Putting these results together, standard analysis of 
£i[m z ] shows, that for sufficient small g and sufficiently 
large 611 a continuous Pi — Cj transition can take place 

at b^c\g), determined by precisely which of the above 
regimes it falls into. Because the driving force (the ex- 
change energy) behind the canting transition vanishes in 
the by ^> 1 limit, Eq.fj], we predict that P-C transition 
must be reentrant, as illustrated in Fig.|l]. 

So far we have restricted our analysis to a spatially 
uniform mean-field treatment. While this is appropriate 
for the commensurate state and the Pc — Cc transition, 
this clearly fails in the incommensurate state, where the 
soliton array breaks translational symmetry. A more de- 
tailed analysis shows that even within mean-field the- 
ory the canting instability is periodically modulated by 
(d x (j) s (x)) 2 , with canting (m z ^ 0) confined to regions be- 
tween solitons, where the twisting of <j)(x) and therefore 
the exchange energy cost is highest. One consequence of 
this is an upward shift of b^ (g) , corresponding to the ex- 
change energy density cost n s ^/ p sz e c i 2 associated with 
deformation of m z (x) localized on solitons. 

It is noteworthy, that, because e c vanishes while p s 
saturates in the d — > limit 0, in principle, for suffi- 
ciently small intcrlayer separation d, the P—C transition 
must always take place. Whether it remains continuous, 
beyond our above mean-field theory analysis is a more 
difficult question. However, estimates of our model pa- 
rameters from recent experiments 0], suggest that cur- 
rently available bilayer devices have g > j c , and there- 
fore should not display the P — C transition (see Fig.Rl) . 



Transcending these semi-microscopic, model-specific 
mean- field considerations, (±m z ) 1i symmetry dictates 
the form of the effective classical Hamiltonian 



H 



d 2 r 



a 2 2 , 4 

— m z + a 4 m z 



(8) 



valid near the finite temperature Pj — Cc transition and 
on scales longer than the soliton lattice spacing nj 1 . 
Given our discussion above, we expect the reduced tem- 
perature a2 ~ T/T c (bu,g) — 1, the compressional elas- 
tic constant c x to vanish as b\\ — > I + (with precise 
form depending on the range of fry), the tilt modulus 
Cy ~ Ps^s/A, with both approaching p s in the dense 
regime, and w a (pseudo)magneto-elastic coupling of the 
soliton lattice phonon u degree of freedom to the local 
electric dipole moment m z ; because a decrease in the 
soliton density (d x u > 0) increases the exchange energy 
of the twisted incommensurate state, we expect w > 0. 

Standard analysis of H in Eq.^ predicts that both 0J4 
and w are relevant couplings for d < 4 and will therefore 
lead to non-meanfield critical behavior sufficiently close 
to the Pi — Ci transition. Preliminary renormalization 
group (RG) analysis in d = 4— e dimensions suggests that 
the magneto-elastic coupling w can either drive the Ising 
transition first order or qualitatively modify it into a new 
scalar compressible Ising model universality class. jl4 15 
If latter scenario survives down to d = 2, critical behav- 
ior of the magnetization m z and the associated magnetic 
susceptibility leads to the predictions of Eqs.|l|,|| for the 
corresponding interlayer charge imbalance q and the dif- 
ferential capacitance C, with gate voltage V playing the 
role of the associated z-directed "magnetic field" . 

We now turn to dynamics. Adapting the imaginary 
time action of the pseudo-ferromagnet |8J to the incom- 
mensurate soliton state, we find 



S 



h/T 



drd 2 r 



ijm z d T u + H[u, m z ] 



(9) 



where the dynamics originates from the Berry's phase 
"pg" term, that encodes the ferromagnetic precessional 
dynamics, 7 = hn s /2i 2 , and oli cx bi*\g) — 611 at T = 0. 

The dynamics of interlayer charge and soliton lat- 
tice fluctuations contained in S, Eq.^| can be probed 
through a linear response of interlayer charge imbalance 
eSn(k, lu) to a time-dependent interlayer voltage V(k, lu), 
applied in a balanced (i.e., keeping vt = 1) capacitive 
geometry. The relevant response function is the dy- 
namic dielectric constant e(k,o>) = ^(4^2) Xzz(k, lu), 
related to interlayer capacitance via a standard rela- 
tion C(k, to) — e(k,ui)A/d, both expressible in terms 
of the pseudo-spin linear susceptibility, x 2Z (k,w) = 
-i J dtd 2 re- iut -^- r ([m z (r,t),m z (0,0)}) . 
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It is straightforward to compute the e(k, to) away from 



the Pi — Ci transition. For the Pj phase we find 



ep(k, uj) = d 



47r£ 2 7 



c yk y 



-co' 



^p(k) 2 



(10a) 



w P (k) - l-\c x k 2 x + cyklf^iplk 2 + a,) 1 / 2 , (10b) 

a result resembling that for the A = B\\ =0 bilayer ||, 
here, with the translational Goldstone mode u(r,t) play- 
ing the role analogous to the staggered U(l) charge Gold- 
stone phason mode (j)(r,t). It is noteworthy that simi- 
lar to recent theoretical predictions |)| and experimental 
findings ||,[l| for interlayer tunneling, at long wavelengths 
the peak in the dielectric response e(k, uo) traces out the 
soliton lattice phonon dispersion wp(k). 
Inside the Ci phase we instead find 



ec(k,w) = d 



Cxk 2 + Cyk, 



\A*Pl) _( w _2(Jsi)i fe )2 + Wc(k )2' 

(11a) 

Wc (k) = 7 _1 (c*^ + Cyk 2 y )V 2 (p* s k 2 + 2|a 2 |) 1 /2, (lib) 

with e(k — > 0, — * 0) and therefore the associated ca- 
pacitance (cf. Eq.||) diverging as the Pj — Cj transition 
is approached from above or below. 

Near the Pj — C\ transition, the divergent correlation 
length \ mz = \J p%jai leads to breakdown of perturba- 
tion theory, and a full RG analysis of this quantum com- 
pressible Ising transition is in principle necessary. On 
general grounds, sufficiently close to b c (g), we expect 
e(k, u>) to display critical scaling 



XK") = k- 2+r <e(k x ./k z y *,cj/k* 



, kyOi 



2 V ) 



(12) 



where v and z x ^ are the correlation length and 
anisotropy exponents, respectively. Preliminary RG 
analysis fl4j| indicates that the upper critical dimension, 
below which the quartic (a^) and magneto-elastic (w) 
nonlinearities become qualitatively important is d uc = 2. 
Hence, we expect mean- field description with z x = 1, 
z u = 2, rj = 0, and v = 1/2 to accurately (up to logarith- 
mic corrections) describe the Pj — Ci transition in these 
2d devices. 

Up to now we have ignored long-range part of the elec- 
trostatic dipolar interaction J, U(k)\m z (k)\ 2 , with kernel 

U(k) = fff(l - e~ kd )/k, and s d = (£^f ^ is the 
dipolar energy per unit of length. As in real bulk mag- 
nets this interaction clearly favors a development of anti- 
aligned dipolar domains at the shortest possible length 
scales and therefore competes with the exchange energy 
■^Pg|Vm z | 2 , minimized by spatially homogeneous m z (r). 
In contrast to 3d systems, where domain size scales as 
\TL % it is easy to show that in our 2d geometry, suffi- 
ciently below the Pi — Ci transition, the domain length 
is £d efe £cxch / £d , (e cxc h. ~ P z s /^m z ) and therefore is a 
length that can in principle be tuned from bilayer thick- 
ness to a macroscopic length. Sufficiently close to the 



transition, the long-range dipolar interaction will always 
become important, and will lead to the canting insta- 
bility in ui z (r) to take place at a finite wavevector kd, 
determined by ^ = 3^ [1 _ (1 + k d d) e - k * d ] . Al- 
though the mean-field predictions of the thermodynamics 
at the transition should remain unchanged, the uniform 
(k x = 0) pseudo-spin susceptibility, determining the ca- 
pacitance C will no longer diverge. For sufficiently large 
domains we expect fluctuations enhancement of C 
near the transition, but with its divergence now cutoff 
by the dipolar domain size, and therefore with its peak 
scaling according to k d ^ 2 v ^ z *. 

We have argued that in bilayer QH ferromagnets, for 
sufficiently strong field and close interlayer separation a 
quantum interlayer charging transition must take place. 
While arguments in favor of the transition and its striking 
properties are quite general, it is unlikely that our mean- 
field estimates of the phase boundary are quantitatively 
trustworthy. Microscopic calculations, exact diagonaliza- 
tion and quantum Monte Carlo studies are necessary to 
accurately determine the details of the phase diagram 
proposed here. We hope that the signatures of the P — C 
transition studied here will stimulate experimental efforts 
to develop bilayer devices, where it can be observed. 
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